Exploring the Low-Z Shore of the Island of Inversion at N = 19 
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The technique of invariant mass spectroscopy has been used to measure, for the first time, the ground state 
energy of neutron-unbound 28 F, determined to be a resonance in the 27 F + /i continuum at 220(50) keV. States in 
28 F were populated by the reactions of a 62 MeV/u 29 Ne beam impinging on a 288 mg/cm 2 beryllium target. The 
measured 28 F ground state energy is in good agreement with USDA/USDB shell model predictions, indicating 
that pf shell intruder configurations play only a small role in the ground state structure of 28 F and establishing 
a low-Z boundary of the island of inversion for N = 19 isotones. 
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A hallmark of the nuclear shell model is its reproduction of 
large energy gaps at nucleon numbers 2, 8, 20, 28, 50, 82, and 
126. Although well established in stable nuclei, these magic 
numbers begin to disappear for nuclei far from stability. For 
example, it has been known for over 30 years that the large 
shell gap at TV = 20 diminishes for neutron-rich nuclei [1—4]. 
The change in shell structure around TV = 20 is now known to 
be a result of the tensor force, which is strongly attractive for 
the spin flip pairs nd 5 / 2 -vd^/ 2 an d strongly repulsive for the 
pairs ?T^5/2-v/7/2 [5-7]. For nuclei in the region of TV ~ 20 
and Z < 13, the reduced TV = 20 gap allows pf shell intruder 
configurations, in the form of multi-particle, multi-hole (np- 
nh or nHco) cross shell excitations, to compete with standard 
sd only configurations if the gain in correlation energy is on 
the same order as the size of the shell gap [8-10]. This has 
led to the establishment of the "island of inversion" — a region 
of nuclei near TV = 20 for which the intruder configuration is 
dominant in the ground state. 

The island of inversion was originally thought to only in- 
clude those nuclei with 10 < Z < 12 and 20 < TV < 22 [11]. 
In more recent years, it has become clear that the island ex- 
tends further, and much experimental effort has been put forth 
to determine its boundaries [12]. On the low-TV (western) 
and high-Z (northern) sides of the island — both in the direc- 
tion of increasing stability — it is generally agreed that ground 
state intruder components fade away for Z > 13 and TV < 18. 
Ground state observables for nuclei lying outside these lim- 
its are well described by sd shell model calculations, such as 
those utilizing the USDA/USDB effective interactions [13]. 
Heading away from stability, the role of intruder configura- 



tions in nuclear ground states becomes less clear. On the high- 
TV (eastern) side, there are strong indications that ground state 
intruder dominance persists in heavier isotopes of Mg, Na, and 
Ne. For example, recent measurements of 2p knockout cross 
sections from 38 Si to 36 Mg have indicated a Oho ground state 
occupation of only 38(8)% in 36 Mg [14], extending the region 
of inversion to at least TV = 24 for the Mg isotopic chain. 

Until now, the low-Z (southern) side of the island has been 
almost completely unexplored. A measurement of bound ex- 
cited states in 27 F, which lies on the island's western border 
at TV = 18, has hinted at pf shell contributions to its excited 
state structure [15], but mass measurements [16] indicate that 
the F ground state is primarily sd shell. For the heavier 
(TV > 19) fluorine isotopes, lying within the island's western 
boundary, no direct experimental information is available. All 
that is known are dripline systematics, e.g. as outlined in [17]. 
The fluorine isotopic chain is also the only area in which the 
intruder structure of Z < 10, TV > 19 nuclei can be examined 
in any detail since the neutron dripline for all lighter elements 
is located at TV < 16 [18-22]. This means that the TV = 19 iso- 
tones of oxygen and below are unbound with respect to three 
or more neutrons, making their study extremely difficult, if 
not impossible. 

In this letter, we report on the first experimental investiga- 
tion of the low-Z border of the island of inversion for TV > 19 
nuclei. This is done via a measurement of the ground state 
of neutron-unbound 28 F, performed with the technique of in- 
variant mass spectroscopy. Our results are then compared to 
binding energy (mass defect) predictions of the USDA/USDB 
shell model. In particular, investigation of the TV = 19 binding 
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energy systematics provides strong evidence that only the iso- 
tones at Z = 10, 11, and 12 are located within the island. In 
addition to mapping the island of inversion, the evolution of 
intruder structure in neutron-rich fluorine isotopes is relevant 
to the abrupt shift in the neutron dripline observed between 
oxygen (ending at N = 16) and fluorine (ending at N > 22) 
[23]. 

The experiment was performed at the National Supercon- 
ducting Cyclotron Laboratory (NSCL) at Michigan State Uni- 
versity, using a primary beam of 48 Ca 20+ accelerated to 
140 MeV/u in the coupled K500-K1200 cyclotrons [24]. The 
48 Ca beam was fragmented in a 1316 mg/cm 2 beryllium pro- 
duction target, and fragmentation products were selected in 
the A 1900 fragment separator [25], with the third and fourth 
segments set to a rigidity of 3.47 Tm and a momentum accep- 
tance of 3.9% to optimize the transmission of 29 Ne fragments. 
The secondary beam was passed through a pair of plastic tim- 
ing scintillators, the first located at the A 1900 focal plane and 
the second 44.3 cm upstream of a secondary reaction target. 
Additionally, the beam was sent through a pair of position 
sensitive Cathode Readout Drift Chambers (CRDCs) and a 
focusing quadrupole triplet. The desired 29 Ne composed ap- 
proximately 2% of the beam and could be fully separated from 
other components using time of flight (ToF) and energy loss 
measurements. The incoming rate of 29 Ne beam particles was 
approximately 70 s , and the median energy of the 29 Ne was 
62 MeV/u. 

The secondary reaction target was beryllium with a thick- 
ness of 288 mg/cm 2 . States in neutron-unbound 28 F were pop- 
ulated by one-proton knockout and decayed by neutron emis- 
sion to 27 F + « with a timescale on the order of 10~ 21 s. 
The decays of the unbound states potentially feed bound ex- 
cited states in 27 F, necessitating the measurement of neutrons, 
charged fragments, and y-rays. The gammas were detected 
using the CAESAR CsI(Na) array [26], which surrounded the 
target and provided an in-beam detection efficiency of ~ 30% 
for 1 MeV gammas. Charged fragments were deflected 43° by 
the Sweeper magnet [27] and passed through a pair of CRDCs, 
an ionization chamber, and thin (0.5 cm) and thick (15 cm) 
plastic scintillators. Neutrons were detected in the Modular 
Neutron Array (MoNA) [28], whose front face was located 
658 cm downstream of the target at 0°. MoNA was partially 
shadowed by the entrance flange of the Sweeper vacuum box, 
resulting in a neutron angular acceptance of roughly ±6.0° in 
the dispersive (x) plane and ±2.5° in the non-dispersive (y) 
plane. 

The charged particle measurements allowed for event-by- 
event isotope selection. Elements were selected using cuts on 
both AE-ToF and AE-E, with AE taken from the ion cham- 
ber signal and E taken from the total charge collected in the 
thick scintillator. For a given element, isotopes were separated 
by correcting their ToF for the various paths taken through 
the Sweeper. This amounts to removing correlations between 
ToF and a variety of other measured parameters, the most 
important being the dispersive position and angle exiting the 
Sweeper (up to fourth order), non-dispersive position exiting 
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FIG. 1 . (Color online) Corrected ToF for fluorine isotopes produced 
from 29 Ne. 
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FIG. 2. (Color online) Simulated resolution and acceptance of the 
experimental setup. Each colored histogram was generated by simu- 
lating a 28 F breakup at the indicated relative energy, then folding in 
detector resolution and acceptance cuts. The shaded curve was gen- 
erated by simulating a 28 F breakup with the relative energy uniformly 
distributed from 0-3 MeV and folding in acceptance and resolution. 
The colored histograms are all normalized to a total area of unity, and 
the shaded curve was arbitrarily scaled to fit within the same panel. 

the magnet, and dispersive position of the incoming beam. A 
plot of the corrected ToF for fluorine elements is shown in 
Fig. 1, with 27 F indicated. 

The decay energy, E c /, of the breakup of unbound states was 
calculated using invariant mass analysis: 



Ed = y mj + ml + 2 (E f E„ - p f p n cos 9) - m f -m„, ( 1 ) 

where mf(m n ), Ef(E n ), and Pf(p n ) refer to the mass, en- 
ergy, and momentum of the charged fragment (neutron), re- 
spectively, and 6 is the opening angle between the two decay 
products. The neutron input to Eq. (1) was calculated from 
ToF and position measurements in MoNA using linear kine- 
matics, while the charged fragment input was reconstructed 
from measurements of the post-Sweeper emittance and the x 
position of the beam on target [29]. 
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FIG. 3. (Color online) Measured decay energy spectrum (including 
smearing from experimental resolution and acceptance) for 27 F + n 
coincidences. The filled squares with error bars are the measured 
data, and the dashed red and dotted blue curves represent the 220 
keV and 810 keV simulation results, respectively. The solid black 
curve is the sum of the two resonances, with the ratio of 220 keV 
resonance to the total area being 28%. The filled orange curve is a 
simulation of a single resonance at 590 keV, and the grey dot-dashed 
curve is the best fit of a single .s-wave (a s = —0.05 fm). 

Resonant states were modeled by a Breit-Wigner lineshape 
with energy dependent width, derived from /^-Matrix theory 
[30]. Free parameters are the central resonance energy Eq, the 
resonance width Fq, and the relative contribution to the overall 
decay spectrum. The orbital angular momentum was fixed at 
£ = 2 assuming the emission of a (W3 /2 neutron. This assump- 
tion may be incorrect if intruder components are significant in 
28 F; however, separate analyses of the data using £ = 1 and 
£ = 3 resonances give results that do not differ significantly 
from the £ = 2 case. We also investigated the possibility of 
£ = decay by modelling the data as an s-wave [31] with the 
scattering length, a s , freely varying. 

Smearing from experimental resolution and acceptance was 
accounted for in a Monte Carlo simulation of the experimen- 
tal setup, including all relevant detector resolutions and accep- 
tance cuts. The simulation included a realistic incoming beam 
profile and used the Goldhaber model [32] to describe the one- 
proton knockout reaction populating 28 F. A demonstration of 
the simulated overall resolution and acceptance functions is 
shown in Fig. 2. To determine optimal fit parameters, large 
Monte Carlo data sets (~ 3 million events) were generated and 
compared to the experimental data using an unbinned maxi- 
mum likelihood technique [33]. 

The measured decay energy of 28 F into 27 F + n (with no un- 
folding of resolution or acceptance) is shown in Fig. 3. Com- 
parison with Fig. 2 reveals that the measured data are strongly 
distorted by resolution and acceptance. In particular, the width 
of the measured data is almost entirely due to experimental 
resolution, and the shape of the data above ~ 0.8 MeV is dom- 
inated by the limited acceptance at higher relative energies. 
No coincident gamma events were recorded in CAESAR, so 
the observed transitions are assumed to feed the ground state 



of F (around 30 CAESAR counts would be expected in the 
case of 100% branching to an exited state). States in 28 F were 
populated in a direct one-proton knockout reaction from 29 Ne, 
so non-resonant contributions to the data are not expected. An 
attempt to fit the data with a single Breit Wigner resonance is 
shown as the shaded orange curve in the figure. In this fit, 
the width was limited to Yq < 1 MeV, beyond which the func- 
tion saturates and lineshapes become indistinguishable. The 
optimal one-resonance fit parameters are Eq = 590 keV and 
To = 1 MeV. The To value is two orders of magnitude larger 
than the single particle prediction, already suggesting that a 
single resonance cannot properly describe the data. Further- 
more, a visual comparison of simulation and data indicates a 
poor fit, with the simulated curve being significantly narrower 
than the data despite its unphysically large Fq value. The best- 
fit of a single s-wave (a s = —0.05 fm) is shown as the grey 
dot-dashed curve and is clearly not consistent with the data. 

The poor fit of a single resonance suggests that multiple 
resonances are present in the data, and a good agreement be- 
tween simulation and data is achieved by modelling with two 
independent Breit Wigner resonances, with the width of each 
resonance fixed at its approximate single-particle value. Given 
this model, the best fit is obtained with the lower resonance at 
220(50) keV (r = 10 keV) and the upper resonance at 810 
keV (ro = 100 keV), as shown in Fig. 3. The relative con- 
tributions of the lower and upper resonances to the total area 
are 28% and 72%, respectively. The stronger population of an 
excited state is somewhat surprising; however, as discussed 
below, contributions from multiple unresolved excited state 
resonances may be present in the data. Thus it is not possible 
to draw any definite conclusions from the relative areas of the 
two resonances. To examine the statistical significance of the 
two-resonance versus one-resonance hypotheses, we calculate 
the likelihood ratio (— 21n[Li /L2] , where L\ and L2 are the re- 
spective likelihoods of the one-resonance and two-resonance 
hypotheses) to be 22.6, providing a moderate level of support 
for the two-resonance model. 

As mentioned, it is possible that more than two resonances 
are present in the data, but the low statistics and limited ex- 
perimental resolution prevent the inclusion of additional res- 
onances into the fit with any sort of certainty regarding their 
location. However, attempts to fit the data with three or more 
resonances demonstrate that the location of the ground state 
peak is insensitive to the presence of multiple resonances. It 
should be noted that the placement of the ground state peak 
at 220 keV is primarily dictated by the shape of the spectrum 
below ~ 300 keV and not by the dip at around 400 keV. In 
particular, modelling with a strong resonance at lower decay 
energies (< 150 keV) results in a sharp peak that is not consis- 
tent with the gradual rise from zero seen in the data. Likewise, 
placing the lowest resonance too high results in a model that 
significantly under-predicts the number of events below ~ 300 
keV. 

Combining the present measurement of the 28 F neutron sep- 
aration energy with the mass measurements of [16], which 
place the 27 F atomic mass excess at 24630(190) keV, we de- 
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FIG. 4. (Color online) Difference between experimental and theoret- 
ical (USDA, USDB) binding energies for N = 19 isotones, 9 < Z < 
17. The error bars on the data points represent experimental errors 
only. The blue dotted, red dashed, and black dash-dotted bands rep- 
resent the respective 170 and 130 keV RMS deviations of USDA and 
USDB interactions. Experimental values, save for Z = 9 which is 
from the present work, are taken from [16] if reported there; other- 
wise they are from the 2003 Atomic Mass Evaluation [34]. 

termine the 28 F binding energy to be 186040(200) keV. As 
presented in [13], it is possible to deduce the presence of 
ground state intruder components in < 20 nuclei by com- 
paring experimental binding energies to sd shell model pre- 
dictions, such as those of the USDA and USDB effective in- 
teractions. For a given nucleus, good agreement between ex- 
periment and USDA/USDB theory indicates a ground state 
configuration that is primarily sd shell. In contrast, a nu- 
cleus with significant ground state intruder components will 
be poorly described by the USDA/USDB shell model. 

Fig. 4 presents a plot of BE exp — BE^ for N — 19 iso- 
tones, 9 < Z < 17, with the fluorine data point taken from 
the present work. The agreement between experiment and 
USDA/USDB predictions is good for the isotones closer to 
stability (Z > 13). For Z = 10-12, the USDA/USDB calcu- 
lations predict significantly lower binding than experiment, as 
expected for these island of inversion nuclei. For 28 F, the good 
agreement between experiment and USDA/USDB is recov- 
ered, providing evidence that intruder configurations are not 
significant in the ground state of 28 F. 

In conclusion, we have determined, for the first time, the 
ground state of 28 F to be a resonance 220(50) keV above the 
ground state of 27 F using the technique of invariant mass spec- 
troscopy. Combined with the mass measurements of [ 16], this 
translates to a 28 F binding energy of 186040(200) keV. Inves- 
tigation of N = 19 binding energy systematics, including the 
present measurement, shows good agreement between exper- 
iment and USDA/USDB predictions for 28 F, in sharp contrast 
to the island of inversion nuclei 29 Ne, 30 Na, and 31 Mg. This 
indicates that pf shell intruder components play only a small 
role in the ground state structure of 28 F, establishing a "south- 
ern shore" of the island of inversion. 

The authors acknowledge the efforts of NSCL operations 
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